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Abstract 

We investigate in this note the dynamics of a one-dimensional KeUer-Segel type model on 
the half-line. On the contrary to the classical configuration, the chemical production term is 
located on the boundary. We prove, under suitable assumptions, the following dichotomy which 
is reminiscent of the two-dimensional Keller-Segel system. Solutions are global if the mass is 
below the critical mass, they blow-up in finite time above the critical mass, and they converge to 
some equilibrium at the critical mass. Entropy techniques are presented which aim at providing 
quantitative convergence results for the subcritical case. This note is completed with a brief 
introduction to a more realistic model (still one-dimensional). 

Resume 

Nous etudions dans cette note la dynamique d'un modele unidimensionnel de type Keller-Segel 
pose sur une demi-droite. Dans le cas present, la production du signal chimique est localisee sur le 
bord, au lieu d'etre repartie a I'interieur du domaine comme dans le cas classique. On demontre, 
sous des hypotheses convenables, la dichotomic suivante qui rappelle le systeme de Keller-Segel en 
dimension deux d'espace. Les solutions sont globales si la masse est sous-critique, elles explosent 
en temps fini si la masse depasse la masse critique. Enfin, les solutions convergent vers un etat 
d'equilibre lorsque la masse est egale a la valeur critique. Des methodes d'entropie sont developpees, 
dans le but d'obtenir des resultats de convergence quantitatifs. Cette note est enrichie d'une breve 
introduction a un modele plus realiste (a nouveau unidimensionnel). 

Version frangaise abregee 

Dans cette note nous allons etudier le coniportement mathematique en dimension un de I'equation 
aux derivees partielles suivante : 

dtn{t,x) — dxxn{t,x) + n{t,Q)dxn{t,x) , i > , a: e (0, +00) , (1) 

avec la condition initiale : n{t = 0, x) = no(a;) > 0. Nous imposons au bord une condition de flux nul : 
dxn{t, 0)+n{t, Oy = 0, de sorte que la masse est conservee au cours du temps (au moins formellement) : 

n{t,x)dx^ / no{x)dx = M. (2) 

2;>0 J x>0 
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Ce modele a ete propose dans [IS] pour decrire synthetiquement la polarisation des cellules de 
levure. Une caracteristique interessante de ([U reside dans le fait que la solution pent devenir non 
bornee en temps fini. Dans cette note nous allons montrer I'alternative suivante : 

Theoreme 1 (Existence globale vs. explosion) Supposons que no{x) est continue sur [0, +00) et 
que riQ G L}^_{{l + x)dx). Si M < 1 alors la solution de ([1]) est globale en temps. Au contraire si M > 1, 
en supposant en outre que uq est decroissante, alors la solution de ^ explose en temps fini. 

Nous annongons egalement les resultats suivants concernant le comportement asymptotique de la 
solution lorsque M <1 : 

Theoreme 2 (Comportement asymptotique) Dans le cas critique M — 1, il existe une famille 
d'etats stationnaires pour JT]) parametree par a > 0. La solution converge (au sens de I'entropie relative 
(I12p ) vers I'equilibre tel que = J^^^xuq^x) dx. 

Dans le cas sous-critique M < 1, la solution decroit vers zero, et converge (au sens de I'entropie 
relative) vers un unique profil auto-similaire. 

Enfin, nous nous interessons a I'etude d'un modele plus realiste, qui prend en compte I'echange entre 
des particules libres, et des particules fixees au bord qui creent le potentiel attractif (concentration 

m ■■ 

( dtn{t, x) = dxxn{t, x) + ^{t)dxn(t, x) , t > , x e (0, +00) , 
\ Ai'(t) =n(i,0)-/i(t), 

avec la condition de flux au bord : dxn{t, 0) + ^{t)n{t, 0) = /i'(t). 

Theoreme 3 Avec les hypotheses des theoremes precedents, et dans le cas sur-critique M > 1, fi{t) 
converge vers Ji ^ M — \ et la densite n{t, x) converge en entropie relative vers h{x) — /lexp {—Jix). 



English Version 

In this note we shall study the mathematical behavior of the following one dimensional partial differ- 
ential equation: 

dtn{t,x) = dxxn{t,x) + n{t,Q)dxn{t,x) , t > , x € (0, +00) , (3) 

together with the initial condition: n{t — Q,x) — no{x) > 0. We impose a zero-flux boundary condition 
for the density n, 

dxn{t,0) +n{t,Of =0. (4) 

Notice that ([4]) and hq G L\_ guarantees nonnegative solutions n(t, x) > Q and mass conservation (at 
least formally): 

/ n{t,x)dx— I nQ{x)dx — M. (5) 

Jx>0 J x>0 

This model has been proposed in [25] to describe basically the polarisation of cells. The interesting 
feature of ([3]) is that the solution may become unbounded in flnite time. Such a behavior is called 
blow-up in finite time. In this note we shall prove the following simple alternative: 



2 



Theorem 1 Assume uq^x) is continuous on [0, +cxi) and uq G -^+((1 + x)dx). If M < 1 the solution 
of ©-(Ul) is global in time. On the contrary if M > 1, assume in addition that no is non increasing, 
then the solution of blows-up infinite time. 

Remark 2 It would be possible to weaken the assumptions on no (a;) (basically J^^q nQ{x) \\ognQ{x)\ dx < 
+00 ) by using strong regularizing effects of the laplacian ( at least in the subcritical case AI < 1) but 
this is beyond the scope of this note. 

Remark 3 Such a critical mass phenomenon (global existence versus blow-up depending on the initial 
mass) has been widely studied for the Keller-Segel (KS) system: dtn(t, x) = dxx'n{t, x)—dx(n(t, x)dxc(t, x)) 
(also known as the Smoluchowski-Poisson system) in two dimensions of space (see and references 
therein). The KS system describes macroscopically a population of diffusive particles which attract each 
other through a diffusive chemical signal (resp. gravitational field), solution of the Poisson equation: 
—Ac{t,x) = n{t,x) with homogeneous Neumann boundary conditions \2(A fTM Wl[ [7^ . On the other 
hand the chemical field in ^ is in fact solution of the Laplace equation with inhomogeneous Neu- 
mann boundary conditions: —dxc{t, 0) = n{t, 0) (production of the signal is located on the boundary) . 
Although the Keller-Segel cannot exhibit blowing-up solutions in one dimension of space, it is indeed 
the case for ([3]) (Theorem\^. As a conclusion, ([3]) appears to have the same ''singularity" as the 
two-dimensional Keller-Segel system. Note that there exist other ways to mimick the two dimensional 
case singular behaviour of KS in one dimension \1(A 

Remark 4 There is a strong connection between the equation under interest here ([3]) and the one- 
dimensional Stefan problem. The later writes indeed U5}j : 

( dtu{t, z) = dzzu{t, z) , t > 0, z € (-oo,s(t)) , 

\ lim,^_oo 9zw(t, = , u{t,s{t))^0, d,u{t,s{t))^~s'{t). ^> 

The temperature is initially nonnegative: u{0,z) = Uo{z) > 0. By performing the following change of 
variables: (j){t,x) = —u{t,s{t) — x), we get the following equation: 

{ dt(j>{t,x) = dxx4>{t,x) - s'{t)dx(l){t,x) , i > , a; e (0, +00) , , . 

\Ymix^+^dxct^{t,x)^Q, Ht,0)^0, dxct^{t,Q) ^ -s' [t) . ^' > 

By differentiating this equation, we recover ^ forn{t,x) = dx(t>it,x). The condition (j){t,0) — turns 
out to be the mass conservation ofn(t,x). 

This connection provides some insights concerning the possible continuation of solutions after blow- 
up ]15\j . This question has raised a lot of interest in the past recent years J 161 \28\. \24\ \14^ . It is postulated 
in J15\l that the one- dimensional Stefan problem is generically non continuable after the blow-up time. 

Using ad-hoc entropy methods (which are to be adapted to the nonlinearity in this problem), we are 
able to investigate long-time behaviour in the critical (M — 1) and the subcritical case (M < 1): this 
is the purpose of Theorems [5] and [51 In short, the results read as follows: there exists a one-parameter 
family of stationary states for the critical mass only (namely decreasing exponentials). In this case 
the conservation of the first momentum enables to select one particular profile among this family. In 
the subcritical case, an appropriate rescaling has to be performed in order to capture the intermediate 
asymptotics. For each mass M < 1 there exists a unique stationary state (with explicit formulation), 
and we prove convergence (in relative entropy) of the rescaled solution towards this profile (namely 
the product of a decreasing exponential and a Gaussian profile). The results are clearly similar to the 
classical Keller-Segel in two dimensions [S], except that the density converges towards a Dirac mass in 
the critical case [5]. 
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1 The critical mass phenomenon 



Blow-up for M > 1. To prove that solutions blow-up in finite time, we show that the first momentum 
of n{t,x) cannot remain positive for all time. This technique was first used by Nagai [H], then by 
many authors in various contexts (see [H [3l [121 HSl [H] for instance) . Other strategies have been used 
to prove the existence of blowing-up solutions (either constructive by Herrero and Velazquez [16j or 
undirect [17j). however up to date this trick is the only way to provide explicit criterion and appears 
to be quite robust to variations around Keller-Segel [H [8] . 

First, the assumption that no is a nonincreasing function guarantees that n{t,-) is also a non- 
increasing function for any time t > due to the maximum principle (notice that the derivative 
v{t,x) = dxn(t,x) satisfies a parabolic type equation, is initially nonpositive, and is nonpositive on 
the boundary due to (IH)). Therefore —dx'n{t, x)/n{t, 0) is a probability density at any time t > 0. We 
deduce from Jensen's inequality the following interpolation estimate: 

-dxn{t,x) V , /■ 2-dxn(t,x) 
X r — da; < / x — dx , 

< 2n{t, 0) / xn{t, x) dx . 

Jx>0 

Secondly introduce the first momentum 3[t) = J^^q xn{t, x) dx. We have for Af > 1: 

^.n(t,0)-AMt,0)<^(l-Af), (8) 

«!<A/^(1-Af). 
dt 

Therefore blow-up of the solution occurs in finite time if A/ > 1. 

Global existence for M < 1. Global existence results for Keller-Segel type systems have been 
initiated by Jager and Luckhaus 19J in the two dimensional case. It relies on a mixture of Gagliardo- 
Nirenberg type and interpolation inequalities. The novelty here is to use a trace- type Sobolev inequality 
(simple in the one-dimensional setting) which is required due to the location of the chemical source on 
the boundary. 

We compute the evolution of the density entropy as following, 



— / n{t^ x) log 'n{t, x) dx = I dtn{t, x) log n{t, x) dx 

x>0 Jx>0 



dt 



{dxn{t, x) + n(t, 0)n{t, x)) ^^^^^^ da; 
x>o n(t,x) 

{dx logn{t, a;))^ n{t, x) dx + n{t, 0)^ . 

x>Q 

The one-dimensional trace inequality we mentioned above writes as following, 

n{t,Q) — — I dxn{t,x)dx^— I {dxlogn{t^x))n{t,x)dx ^ 

n{t,Of<M {dxlogn{t,x)fn{t,x)dx. (9) 



Therefore we deduce that 
d_ 
di 



x>0 



n{t,x)logn{t,x)dx < {M - 1) I {dxlogn{t,x)f n{t,x) dx , (10) 



a;>0 Jx>0 
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hence the entropy is nonincreasing when the mass is smaUer than 1. Observe that equahty holds in 
the trace inequaUty ^ if \ogn{t,x) is constant w.r.t. x: there exists a(t) > such that n{t,x) = 
Ma(t) exp{—a{t)x). In fact the boundary condition ^ imphes M = 1, which is the only configuration 
where a stationary state can exist (see Section [2]) ■ 

A major step towards a complete existence theory of ([3]) in the subcritical is to ensure that the 
boundary value n{t, 0) makes perfect sense. This is a consequence of an Aubin-Lions type argument 
[T], which is straightforward in this over-simplified context. We ask for continuity w.r.t. x of the 
density n(t, x): 

{n{t,x)^n{t,y)f < (^J\{t,z)dz^ (^£''(9, logn(i, z) dz^ . 

The bound (fTO|) together with the control of moments guarantee that J^^gidx ^ogn{t, x))'^n{t, x) dx is 
finite almost every time. Therefore n{t, •) is continuous almost every time. 

To conclude this Section, let us mention that it is now classical to prove suitable regularizing effects 
acting on ([3]) in the subcritical case M < 1. Indeed an a priori estimate (fTO|) on /^^^ n{t, x) (log n{t, x))^ dx 
yields the boundedness of all L''— norms {1 < p < +00) [TOl [71 llOj. 

2 Long-time behaviour: convergence in relative entropy for 
the critical and the subcritical cases 

The critical case. Equilibrium configurations for the cell density are only possible when the mass 
is critical: M — 1 (as it is for the two-dimensional Keller-Segel problem). In this case, a straigtforward 
computation leads to the one-parameter family: 

ha{x) = ae-""" , a>0. (11) 

On the other hand, notice that the first momentum of the cell density is conserved This prescribes 
a unique choice for a: a^^ ~ J(0) . 

Theorem 5 Assume nQ{x) being as in Theorem\^ and the mass being critical: M — 1. As time goes 
to infinity, the cell density converges (in relative entropy) towards ha{x). 



The convergence proof is based on evaluating the time evolution of the relative entropy, defined as 
follows: 

The precise description of the equality cases for inequality ([9]) enables to perform accurate estimates. 
A direct computation yields the following estimate: 

-^H(t)=-/" {dx\ogn{t,x)+n{t,{))f n{t,x)dx. (13) 
dt Jx>o 

We refer to [5] for more details. 



Self-similar decay in the subcritical case. In the sub-critical case M < 1 one expects the density 
n(t, x) to decay self-similarly. For this purpose the density is appropriately rescaled: 

n(t,x) = — ^ u(\ogy/l + 2t, ) . 
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The new density u(T,y) satisfies: 



dru{T, y) = dyyu[T, y) + dy {yu{T, y)) + u{t, 0)dyu{T, y) , (14) 

and no-flux boundary conditions: dyu{T, 0) + u(r, 0)^ = 0. The additionnal left-sided drift contributes 
to confine the mass in the new frame (r, y) . The stationary equihbrium in this new setting can be 
computed expHcitely. The expected self-similar profile writes: ga{y) — aexp (—ay — 2/^/2), where a 
is given by the relation P{a) = M, P being an increasing function defined as following: 

"^^"^ ^ X>o " ^ j ' I lim«^.oo P(a) ^ 1 ■ 

Theorem 6 Assume no{x) being as in Theorem]^ and the mass being subcritical: M < 1. As time 
goes to infinity, the first momentum J(t) of the density converges to a{l — M) and the cell density 
converges (in relative entropy) towards ga{y). 



H(r)^/ ^logr^U.(.)d.. (15) 



The proof of this Theorem relies again on the time evolution of the relative entropy 

uir, y) / u{t, y] 

'y>Q 9a{y) \ 9a{y) 

More precisely we have: 

A /H(r) + ^ (J(r) - a(l - M)f 

1 / d \ ^ 

^M(T,y)(aj^logw(T,?/) + y + u(T,0))^ dy--^^-— ^ ^— J(r)j . (16) 

3 Analysis of a coupled ODE/PDE model 

We investigate in this section a variant of ^ which is more relevant for modelling purposes [25] . In 
this new setting, the chemical is supplied by a quantity /i(i) which evolves by exchanging particles at 
the boundary x = Q: 

dtn{t,x) ^ dxxn{t,x) + ^{t)dxn{t,x) , t > , x € {0, +oo) , , , 

^l'{t)=n{t,o)-^,{t), ^''> 

together with the initial conditions: n{t = 0,2:) = no{x) > and /i(i = 0) = fiQ. The conservation of 
the total mass of particles: 

n{t,x)dx + n{t) ^ M , (18) 



'x>0 

yields the following boundary condition for the cell density: 

dxn{t, 0) + iJ.{t)n{t, 0) = ii'{t) . 

Long-time convergence in the case M > 1. We denote by m{t) the mass of the cell density 

nit, x): 

m{t) = I n{t,x)<lx. (19) 
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(notice m'{t) + iJ.'{t) — due to the conservation of mass). Introduce the relative entropy: 



Jx 



>o ra{t)h{x) 



n{t, x) 



log 




where the expected profile h is given by: 



h{x) = /i exp {—^x) , /I = A/ — 1 . 



Theorem 7 As time goes to infinity, the mass m{t) of the cell density converges to 1 and the cell 
density converges (in relative entropy) towards h{x). 

The proof of this Theorem relies again on the time evolution of the relative entropy. This is strongly 
inspired from the previous computation, but takes into consideration the non-conservation of mass for 
the cell density and the dynamics of n{t) [9]. 
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